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We present exact solutions of the gravitational field equations in the generalized Randall-Sundrum 
model for an anisotropic brane with Bianchi type I geometry, with a generalized Chaplygin gas as 
matter source. The generalized Chaplygin gas, which interpolates between a high density relativistic 
era and a non-relativistic matter phase, is a popular dark energy candidate. For a Bianchi type 
I space-time brane filled with a cosmological fluid obeying the generalized Chaplygin equation of 
state the general solution of the gravitational field equations can be expressed in an exact parametric 
form, with the comoving volume taken as parameter. In the limiting cases of a stiff cosmological 
fluid, with pressure equal to the energy density, and for a pressureless fluid, the solution of the field 
equations can be expressed in an exact analytical form. The evolution of the scalar field associated 
to the Chaplygin fluid is also considered and the corresponding potential is obtained. The behavior 
of the observationally important parameters like shear, anisotropy and deceleration parameter is 
considered in detail. 

PACS numbers: 04.50.+h, 04.20.Jb, 04.20.Cv, 95.35.+d 



I. INTRODUCTION 

The observations of high redshift supernovae Q] and the Boomerang/Maxima data Q, showing that the location of 
the first acoustic peak in the power spectrum of the microwave background radiation is consistent with the inflationary 
prediction Q — 1, have provided compelling evidence for a net equation of state of the cosmic fluid lying in the range 
— I < w = p/p < —1/3. To explain these observations, two dark components are invoked: the pressureless cold dark 
matter (CDM) and the dark energy (DE) with negative pressure. CDM contributes Sl m ~ 0.3, and is mainly motivated 
by the theoretical interpretation of the galactic rotation curves and large scale structure formation. DE is assumed to 
provide Qde ~ 0.7 and is responsible for the acceleration of the distant type la supernovae. There are a huge number 
of candidates for DE in the literature (for recent reviews see Q and 0). However, neither CDM nor DE has direct 
laboratory observational or experimental evidence for its existence. It would be important if a unified dark matter 
-dark energy scenario could be found, in which these two components are different manifestations of a single fluid 0. 

A candidate for such an unification is the so-called generalized Chaplygin gas, which is an exotic fluid with the 
equation of state p = —B/p a , where B and a are two parameters to be determined. It was initially suggested by 
Kamenshchik, Moschella and Pasquier [|| with a = 1, and then generalized by Bento, Bertolami and Sen 7] for the 
case a 1. The Chaplygin gas also appears in the stabilization of branes in Schwarzschild-AdS black hole bulks as 
a critical theory at the horizon p| and in the stringy analysis of black holes in three dimensions . The Chaplygin 
equation of state can also be derived from Born-Infeld type Lagrangians Q, [H|- This simple and elegant model 
smoothly interpolates between a non-relativistic matter phase (p = 0) and a negative-pressure dark energy dominated 
phase. 

The cosmological implications of the Chaplygin gas model have been intensively investigated in the recent literature 

m 

The Chaplygin gas cosmological model has been constrained by using different cosmological observations, like type la 
supernovae |l3j |. |14|. the CMB anisotropy measurements |l5j . gravitational lensing surveys fl6| . the age measurement 
of high redshift objects ^3 an d the X-ray gas mass fraction of clusters ^8|. The obtained results are somewhat 
controversial, with some of them claiming good agreement between the data and the Chaplygin gas model, while the 
rest ruling it as a feasible candidate for dark matter. In particular, the standard Chaplygin gas model with a = 1 is 
ruled out by the data at a 99% level jig ■ 

The embedding approach to gravity, according to which our four-dimensional space-time is a three-brane, embedded 
in a five-dimensional space-time (the bulk), has a long history, beginning at least with the classic book by Eisenhart 
[T^|. More recently, this idea was considered by Akama [2(j and Rubakov and Shaposhnikov pHl (a n exhaustive list 
of references on the early articles published in this field can be found in Bandos and Kummer |22 ] and Pavsic and 



'Electronic address: mkmak@vtc.edu. hk 
tElectronic address: harko@hkucc.hku.hk 



2 



Tapia |23|. respectively). Recently, due to the proposal by Randall and Sundrum [24| the idea of the embedding of 
our Universe in a higher dimensional space had attracted again a considerable interest. According to the brane- world 
scenario, the physical fields (electromagnetic, Yang-Mills etc.) in our four-dimensional Universe are confined to the 
three brane. These fields are assumed to arise as fluctuations of branes in string theories. Only gravity can freely 
propagate in both the brane and bulk space-times, with the gravitational self-couplings not significantly modified. 
This model originated from the study of a single 3-brane embedded in five dimensions, with the 5D metric given 
by ds 2 = e"^ v ^r\ ilv dx >L dx v + dy 2 , which, due to the appearance of the warp factor, could produce a large hierarchy 
between the scale of particle physics and gravity. Even if the fifth dimension is uncompactified, standard AD gravity is 
reproduced on the brane. Hence this model allows the presence of large, or even infinite non-compact extra dimensions. 
Our brane is identified to a domain wall in a 5-dimensional anti-de Sitter space-time. For a review of dynamics and 
geometry of brane Universes see 25]. 

Due to the correction terms coming from the extra dimensions, significant deviations from the Einstein theory occur 
in brane world models at very high energies 26] . Gravity is largely modified at the electro- weak scale 1 TeV. The 
cosmological implications of the brane world theories have been extensively investigated in the physical literature • 
Gravitational collapse can al so p roduce high energies, with the five dimensional effects playing an important role in 
the formation of black holes [23 • 

The possibility that the non-local effects, generated by the free gravitational field of the bulk in a conformally 
symmetric brane, may provide an explanation for the dynamics of the neutral hydrogen clouds at large distances 
from the galactic center, and hence provide a purely gravitational explanation for the existence of dark matter, was 
suggested in j2(| . 

The behavior and the properties of the anisotropic Bianchi type I brane-world for different type of matter sources, 
including inflationary scalar fields, have been intensively investigated in the framework of the brane world cosmologies 
|30|. The study of anisotropic homogeneous brane world cosmological models have shown an important difference 
between these models and standard four-dimensional general relativity, namely, that brane Universes are born in an 
isotropic state. For Bianchi typ e I and V geometries this type of behavior has been found both by exactly solving the 
gravitational field equations (3J , or from the qualitative analysis of the model ]Q . A general analysis of the anisotropy 
in spatially homogeneous brane world cosmological models has been performed by Coley [33j, who has shown that 
the initial singularity is isotropic, and hence the initial conditions problem is naturally solved. Consequently, close to 
the initial singularity, these models do not exhibit Mixmaster or chaotic-like behavior j34[. Based on the results of 
the study of homogeneous anisotropic cosmological models Coley conjectured that the isotropic singularity could be a 
general feature of brane cosmology. The analytic study of the values of the anisotropy in the limit of small times shows 
that for all Bianchi type space-times filled with a non-zero pressure cosmological fluid, obeying a barotropic equation 
of state, the initial singular state on the brane is isotropic 35] . This result is obtained by assuming that in the limit 
of small times the asymptotic behavior of the scale factors is of Kasner-type. For brane worlds filled with dust or with 
a scalar field, the initial value of the anisotropy coincides in both brane world and standard four-dimensional general 
relativistic cosmologies. 

It is the purpose of this paper to investigate some classes of exact solutions of the gravitational field equations in 
the brane world model for the anisotropic Bianchi type I geometry, by assuming that the matter source on the brane 
consist of a generalized Chaplygin gas. In this case the general solution of the field equations can be obtained in an 
exact parametric form, for arbitrary values of the parameters characterizing the equation of state of the matter. The 
inclusion of the quadratic terms in the energy-momentum tensor of the perfect cosmological fluids leads to major 
changes in the early dynamics of the anisotropic Universe, as compared to the standard general relativistic case. We 
also consider the brane dynamics of the scalar field associated to the Chaplygin case and obtain the corresponding 
self-interaction potential. The behavior of the observationally important physical quantities like anisotropy, shear and 
deceleration parameters is considered in detail for all these models. 

The present paper is organized as follows. The gravitational field equations for a Chaplygin gas on a brane with 
an anisotropic Bianchi type I geometry are written down in Section II. In Section III we present the general solution 
of the field equations. In Section III we discuss and conclude our results. 

II. GEOMETRY, MATTER AND FIELD EQUATIONS ON THE BRANE 

In the 5D space-time the brane-world is located as Y(X I ) = 0, where X 1 . I — 0,1,2,3,4 are 5-dimensional 
coordinates. The effective action in five dimensions is |36| 

S = J d 5 X^g~ 5 (J^Rs - A 5 ) + J d A x^~g (JLK* - A + L mattcr ) , (1) 
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where fc| = 8irG5 is the 5-dimensional gravitational coupling constant, A5 is the vacuum energy in the bulk, A is the 
vacuum energy on the brane and where x^, /j, = 0, 1, 2, 3 are the induced 4-dimensional brane world coordinates. R$ 
is the 5D intrinsic curvature in the bulk and K is the intrinsic curvature on either side of the brane. 

On the 5-dimcnsional space-time (the bulk), with the negative vacuum energy A5 and brane energy-momentum as 
source of the gravitational field, the Einstein field equations are given by 

G I j = ktT u , T I j = -A 5 g I j + d(Y)[~Xg I j + T}J tt ^], (2) 

In this space-time a brane is a fixed point of the symmetry. In the following capital Latin indices run in the range 
0, 4 while Greek indices take the values 0, 3. 

Assuming a metric of the form ds 2 — {njnj + gujdx 1 dx J , with njdx 1 = d\ the unit normal to the x = const, 
hypersurfaces and gjj the induced metric on \ — const, hypersurfaces, the effective four-dimensional gravitational 
equations on the brane (which are the consequence of the Gauss-Codazzi equations) take the form 26] : 

Gfj, v — —Ag^ + k^T^v + k^S^y — E^ v , (3) 

where is the local quadratic energy-momentum correction 

V = ^TT^ - \T» a T„ a + ± 9fiv (3T^T Q/3 - T 2 ) , (4) 

and is the nonlocal effect from the bulk free gravitational filed, transmitted projection of the bulk Weyl tensor 
Ciajb 

E IJ = C IAJB n A n B , E u -> E^STj as X -> 0. (5) 
The four-dimensional cosmological constant, A, and the coupling constant, are given by 

, 4 ( K k 2 \ 2 \ , 2 k^X 

The Einstein equation in the bulk, Codazzi equation, also implies the conservation of the energy momentum tensor 
of the matter on the brane, 

A, V = 0- (7) 

Moreover, the contracted Bianchi identities on the brane imply that the projected Weyl tensor should obey the 
constraint 

D V EJ = k\D v S^. (8) 

Finally, the equations ©, and (jHJ give the complete set field equations for the brane gravitational field. 
For any matter fields (scalar field, perfect fluids, kinetic gases, dissipative fluids etc.) the general form of the brane 
energy-momentum tensor can be covariantly given as |25| 

Tfu, = pu^u v + phjj, v + TTfu/ + 2q (M u u) . (9) 

The decomposition is irreducible for any chosen 4- velocity u^. Here p and p are the energy density and isotropic 
pressure, and h^ v — g^ v + u^u v projects orthogonal to u M . 

The energy flux obeys the relation = <7< M >, and the anisotropic stress obeys 7r M „ = 7r <Ai;y> , where angular 
brackets denote the projected, symmetric and trace-free part: 



V<^> = h/V v , W <liV> 



W aP . (10) 



The symmetry pro per ties of E^ v imply that in general we can decompose it irreducibly with respect to a chosen 
4- velocity field m m as |25| 



Env 



(11) 
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where k = fcs/fei, U is a scalar, Q M a spatial vector and V^v a spatial, symmetric and trace-free tensor. For 
homogeneous models = and V^ u — |3^ |. Hence the only non-zero contribution from the 5-dimensional 
Weyl tensor from the bulk is given by the scalar or "dark radiation" term U. 

We consider that in Eq. @ the heat transfer is zero, = and we also neglect the effect of the possible anisotropic 
stresses (viscosity of the cosmological fluid, magnetic fields etc.) on the dynamics of the Universe. This allows to 
assume that the cosmological fluid on the brane is perfect, without dissipative effects. Then the matter correction 
terms are given by 

S» v = -^p 2 u ll u u + —p(p + 2p)h iiv . (12) 

The line element of a Bianchi type 1 space-time, which generalizes the flat Friedmann- Robertson- Walker metric to 
the anisotropic case, is given by 

ds 2 = -dt 2 + a\ (t) dx 2 + a\ (t) dy 2 + a\ (t) dz 2 . (13) 

We define the following variable: 

3 

V = Y[ ai , (14) 

i=l 
1 3 

H=-J2H t , (15) 
6 i=i 

Hi = ^,i= 1,2,3, (16) 
at 

AiT i = iTi -H,i = 1,2, 3. (17) 

In equations H14fl - I|17|) V is the volume scale factor, Hi,i = 1,2,3 are the directional Hubble parameters and H is 
the mean Hubble parameter. From Eqs. I|14f> and (|16fl we obtain the important relation H — V/3V. 

The physical quantities of observational interest in cosmology are the expansion scalar 9, the mean anisotropy 
parameter A and the shear scalar c 2 , and are defined according to 

9 = 3H, (18) 

(19) 



i— 1 x / 



a 2 = \a ik a* = \ ^ H 2 3H 2 = \ah 2 . (20) 



;=i 



An important observational quantity is the deceleration parameter 

The sign of the deceleration parameter indicates whether the model inflates or not. The positive sign of q corresponds 
to "standard" decelerating models whereas the negative sign indicates inflation. 

We assume that the isotropic pressure p of the cosmological fluid on the brane obeys a modified Chaplygin gas 
equation of state [33 , 

P = 1P-^, (22) 
P 
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where < 7 < 1 and < a < 1. £> is a positive constant. 

When 7=1/3 and the comoving volume of the Universe is small (p — > 00), this equation of state corresponds to a 
radiation dominated era. The main reason for the use of the modified Chaplygin gas equation of state is that it can 
also describe the very high density (7 = 1) and the radiation dominated phase in the evolution of the Universe. When 
the density is small, p — > 0, the equation of state corresponds to a cosmological fluid with negative pressure (the dark 
energy). In this limit, the modified Chaplygin equation of state corresponds to a mixture of cold dark matter and 
dark energy Q. At all stages it shows a mixture of dark energy, dark matter and radiation. For p = (B/ r y) 1 ^ a+1 ^ the 
matter content is pure dust with p = 0. In the following we neglect the effect of baryonic (luminous) matter on the 
dynamics of the Universe, which, in a rigorous approach, should be added separately to the equation of state of the 
cosmological matter. But since the baryonic matter represents only 4 — 5% from the total amount of energy-matter 
in the Universe [3j, its effect on the large scale cosmological dynamics can be ignored. 

The speed of sound v 2 s — dp/ dp in the Chaplygin gas is given by 



7 (1 



ap 

a) . 

P 



(23) 



Using the variables (|14|) - (|17fl . the Einstein gravitational field equations, the Bianchi identity and the evolution 
equation for the nonlocal dark radiation take the form, respectively, 



3H + J2 H ? = A 



&4 



Vdt K ' 2 
^3H{p + p) = 0, 
U + AHU = 0, 



(3 7 + l)p- 



35" 
~ ~P". 
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~ 2A 
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~ 2A 


IP ~ 



B 



2U 



6U 



1.2.3, 



(24) 

(25) 

(26) 
(27) 



Since we have assumed a homogeneous geometry, the effective nonlocal energy density, V^, and Q M vanish on the 
brane, V^ v = and <2 M = 0, respectively. 



III. GENERAL SOLUTION OF THE FIELD EQUATIONS FOR THE BIANCHI TYPE I BRANE 

By using the equation of state of the cosmological matter given by Eq. ill'L'l) , the energy conservation equation Eq. 
(|2*S|l can be integrated exactly, thus giving the time evolution of the energy density of the modified Chaplygin gas 
filled Bianchi type I brane as 



B 



C 



1 + 7 y(l+7)(l+a) 



where C is an arbitrary constant of integration. 
Integrating Eq. (|T7|) yields 



where Uq is an arbitrary constant of integration. 
By adding Eas. (|25|l we find 

Vdt y ' 3V 2 



B ' 


&4 


2 B ' 




~ 2A 


IP r 

P . 



2U 

fc|A : 



(28) 



(29) 



(30) 



or, by using the explicit dependence of the energy density on the comoving volume V and defining a new function 
V = F (V), we can rewrite Eq. 13U|) in the form 



3P 

V = F(V) = 3AV--+vUy-i) 



B 



C 



2A I ' 



B 



C 



l + 7 y(l+<y)(l+a) 



1 + 7 y(i+7)(i+") 
B B 



B 



B 



C 



C 



1 + 7 v ( - 1+ ^ (1+a '> 
^ 1 6U 



1+7 y(l+7)(l+a) 



Vs. 



(31) 
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The general solution of Eq. I|31|) is given by 

t-t 



where C\ and to are arbitrary constant of integration. 
Subtracting Eq. (|30fl from Eqs. (|25[) we obtain 



dV 



2 J F (V) dV + Ci 



(32) 



H i= H+^,i = 1,2,3, 



(33) 



with Ki,i = 1,2,3 constants of integration satisfying the consistency condition %i = 0. From Eq. 133|l . it is 

easy to show that 



Hf = 3H 2 



K 2 
V 2 



(34) 



Therefore for a Bianchi type I induced brane geometry in the presence of the modified Chaplygin gas the general 
solution of the gravitational field equation can be expressed in the following exact parametric form, with V > taken 
as parameter: 



2fF(V) dV + d 

V : 



(35) 



a l0 V3 exp 



K, 



dV 



V^2 J F (V) dV + Ci 



1,2,3, 



(36) 



.4 



3K 2 



2 / F (V) dV + Ci ' 
K 2 



a 2 = 



2V 2 ' 



(37) 



(38) 



3VF (V) 



2 J F (V) dV + d 



(39) 



where a^o, i = 1,2, 3, are arbitrary constant of integration. In order for the general solution of the field equations given 
by Eqs. H35(l - (|39|l be defined for all V > 0, it is necessary to chose the constants a and B, describing the thermodynamic 
properties of the Chaplygin gas and the integration constant C so that the condition 2 J F (V) dV + C\ > holds 
for all V > 0. Taking into account the definition of V, given by Eq. (|14f> . and Eqs. ( I36|) . and since the integration 
constants Ki, i = 1,2,3 satisfy the condition $^ i=1 K% — 0, it follows that the integration constants a^o, i — 1,2,3, 
have to satisfy the condition fL_i a io = !• 

By introducing two new variables r and v by means of the transformations r = \f2>k±B 1 / 2 ( 1+a H and V = Vqv, 
with Vq =constant, by normalizing the values of the arbitrary integration constants so that C/BVq 



(l+ 7 )(l+a) 



4U V ; yi /kiB^^+^X = 1, Ci/3y 2 fc|S 1 /( 1+Q ) = 1, denoting flV(i+«)/A = b, A = 2A/3fc|S 1 /( 1+Q ) and 

f(v) = 



1 + 7 u(l+7)(l+a) ' 



it follows that Eq. (J32J) can be written in the form 



To 



Jl{ x ov 



dv 



bv 



IP 



'} 



(40) 



(41) 



dv + 1 
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where r = VEkiB^^+^to, 

With the use of Eqs. I|28[) . (|33[l and (|35[) . from Eq. i|24[l it follows that the arbitrary integration constants Ki, 
i = 1,2,3 and C\ must satisfy the consistency condition 



r2 = 2Ci 



3 ' 



(42) 



Scalar fields are supposed to play a fundamental role in the evolution of the early universe. The Chaplygin gas 
model can be also described from a field theoretical point of view by introducing a scalar field <j> and a self interacting 
potential U (</)), with the Lagrangian @, 0> E3> Il3> 



L* = \<? - U{<t>). 

The energy density and the pressure associated to the scalar field <fi are given by 



(43) 



+ U (</>) = p, 



(44) 



and 



respectively. 

In view of Eas.l|i3)l. (IHfll and we have 



2 = (l + 7) 



B 



C 



l + 7 y(i+ 7 )(i+o) 



- B 



B 



and 



17(4) = 



(1-7) 



B 



C 



1 + 7 v (1+ ^ (1+a '> 



— B 

+ 2 



B 



(45) 



C 



1 + 7 y(i+7)(i+a) 



C 



l +7 y(l+ 7 )(l+a) 



(46) 



(47) 



respectively. 

By integrating Eq. (|46|l . we obtain 



\ 



(1 + 7) 


B 1 C 


1+Q 


B , C 


1 + c 


l + 7 1" y(l+ 7 )(l+a) 


l+ 7 "T" y(i+T)(i+a) 






2 J F (V) dV + 







dV, 



(48) 



where </> is an arbitrary constant of integration. 



IV. DISCUSSIONS AND FINAL REMARKS 



In the present paper we have considered Chaplygin gas filled Bianchi type I brane worlds in the framework of the 
Randall-Sundrum scenario [24|. The Randall-Sundrum mechanism was originally motivated as a possible mechanism 
for evading Kaluza-Klein compactification by localizing gravity in the presence of an uncompactified extra dimension. 
In this model, the 5-dimensional bulk space-time is assumed to be vacuum except for the presence of a cosmological 
constant. Matter fields on the brane are regarded as responsible for the dynamics of the brane. The quadratic terms 
in the energy-momentum tensor, due to the form of the Gauss-Codazzi equations, may be important in the very 
early stages of the Universe, leading to major changes in the dynamics of the Universe. The general solution of the 
gravitational field equations on the brane can be obtained for both discussed brane geometries in an exact parametric 
form, with the comoving volume element V taken as parameter. 

By obtaining the general solution of the field equations we have assumed the presence of a non-zero effective 
cosmological constant A on the brane. However, the most attractive feature of the Chaplygin gas is that it could 
explain the main observational properties of the Universe without appealing to an effective cosmological constant. To 
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see if this is indeed the case in the anisotropic brane world models, in the following we will assume that on the brane 
A = 0. This constraint imposes the following relation between the vacuum energies in the bulk and on the brane: 

A 5 = -^A 2 (49) 
o 

Therefore for Ao = and for a fixed 7 and a, the dynamics of the anisotropic brane Universes is controlled by the 
parameter b = _B 1 /( 1 + Q ) / A, the ratio of the Chaplygin gas equation of state parameter and the vacuum energy on the 
brane. 

From the equation of state of the Chaplygin gas with 7 = it follows that for the critical values p c and p c of the 
pressure and density the parameter w c = p c / ' p c is given by w c — —B/p" +1 . Evaluating this relation at the present 
time gives B = — w c o/d" +1 . The Chaplygin gas behaves like a cosmological constant for w c q = —1, which fixes the 
constant B as B = p^ 1 = (3H% /8irG) a+ \ where H = 3.24 x 10~ 18 /i s~\ 0.5 < h < 1, is the Hubble constant 0, 
0. The definition of the dimensionless time parameter r = \/3k4B 1 / 2 ( a+1 H then gives r = 3Hot, or t w lO 17 /!" 1 x r 
s" 1 . 

Generally, the obtained solution of the gravitational field equations describes an expanding brane Universe. The 
variation of the expansion scalar 9 as a function of the dimensionless time parameter r is represented in Fig. 1. 
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FIG. 1: Expansion scalar 6 of the Chaplygin gas filled Bianchi type I brane universe as a function of the dimensionless time 
parameter r for 7 = 1/3, a — 0.5 and b — 0.01 (solid curve), b = 0.02 (dotted curve), b — 0.03 (dashed curve) and b = 0.04 
(long dashed curve). We have assumed a vanishing effective cosmological constant on the brane, Ao = 0. The relation between 
t and the cosmic time t is t ~ 10 17 h~ 1 x r s -1 , 0.5 < h < 1. 

For a Chaplygin type cosmological fluid confined on the brane, the early cosmological evolution in anisotropic 
Bianchi type I geometry is fundamentally changed by the inclusion of the terms proportional to the square of the 
energy density and energy density-pressure products. The time variation of the mean anisotropy parameter of the 
Bianchi type I space-time is represented, for 7 — 1/3, a — 0.5 and for different values of the parameter b, in Fig. 
2. At high densities the brane Universe starts its evolution from an isotropic state, with A(to) = 0. The anisotropy 
increases and reaches a maximum value after a finite time interval t max , and for t > i max , the mean anisotropy is a 
monotonically decreasing function, tending to zero in the large time limit. This behavior is in sharp contrast to the 
usual general relativistic evolution (in this case the quadratic contribution to the energy-momentum tensor vanishes), 
in which the Universe is born in a state of maximum anisotropy. 

The time variation of the shear parameter a 2 is represented, as a function of the dimensionless time r and for 
different values of b, in Fig. 3. 

The shear is a monotonically decreasing function of the cosmological time, and in the limit of large t, t — > 00, 
corresponding to the isotropic limit, a 2 tends to zero. In the limit of small times, since a 2 oc V~ 2 , the shear has a 
singular behavior, tending to infinity, lim t ^o a 2 — > 00. 

In the initial stage the evolution of the Chaplygin gas filled Bianchi type I brane universe is non-inflationary, but 
in the large time limit the brane Universe ends in an accelerating stage. In Fig. 4 we present the dynamics of the 
deceleration parameter for different values of b and for 7 = 1/3 and a = 0.5. 
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FIG. 2: Anisotropy parameter A of the Chaplygin gas filled Bianchi type I brane universe as a function of the dimensionless 
time parameter r for 7 = 1/3, a = 0.5 and 6 = 0.01 (solid curve), b — 0.011 (dotted curve), b — 0.012 (dashed curve) and 
b — 0.013 (long dashed curve). We have assumed a vanishing effective cosmological constant on the brane, Ao = 0. The relation 
between r and the cosmic time t is t « 10 17 /i _1 x r s _1 , 0.5 < h < 1. 




FIG. 3: Shear scalar <r 2 of the Chaplygin gas filled Bianchi type I brane universe as a function of the dimensionless time 
parameter r for 7 = 1/3, a = 0.5 and b — 0.01 (solid curve), b = 0.02 (dotted curve), b — 0.03 (dashed curve) and b = 0.04 
(long dashed curve). We have assumed a vanishing effective cosmological constant on the brane, Ao = 0. The relation between 
r and the cosmic time t is t ~ 10 17 h~ x x r s _1 , 0.5 < /i < 1. 



The time variation of the scalar held associated with the Chaplygin gas is represented for fixed a and 7 and different 
values of the parameter b in Fig. 5. The scalar field is a monotonically decreasing function of time. 



The dependence of the potential U{4>) on the scalar field </> is represented in Fig. 6. 



In the limit of high densities, the equation of state of the Chaplygin gas reduces to the linear barotropic equation 
of state of the relativistic matter, p = 7/9. In this case the density dependence on the comoving volume is given by 
p = p a /V 7+1 , pq =constant, and the equation describing the time variation of the comoving volume reduces to 

t-tp= [—= dV (50) 



^SklpoV^ + 3 ^p 2 V-^ + ±§£v 2 /3 + Ci 



In the limit of very high matter densities, which can be described by the strong equation of state for the hot nucleonic 
gas, it is assumed that p = p. It is believed that matter actually behaves in this manner at densities above about ten 
times the nuclear density, that is, at densities greater than 10 17 g/cm 3 and at temperatures T = (pc 2 7 'urad) 1 ^ > 10 13 
K, where a ra( i is the radiation constant [33]. Therefore in this case 7 = 1. In the limit of small V, the term V~ 2 will 
dominate in the integral in Eq. ||SUJ). By maintaining only this term, the integration can be done easily and the time 
dependence of the relevant geometrical and physical parameters can be obtained as 



1/4 



V = ( ™ ) t V2, (51) 
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FIG. 4: Deceleration parameter q of the Chaplygin gas filled Bianchi type I brane universe as a function of the dimensionless 
time parameter r for 7 = 1/3, a = 0.5 and b — 0.01 (solid curve), b = 0.015 (dotted curve), b — 0.02 (dashed curve) and 
b — 0.025 (long dashed curve). We have assumed a vanishing effective cosmological constant on the brane, Ao = 0. The relation 
between r and the cosmic time t is t » 10 17 ft _1 x r s _1 , 0.5 < h < 1. 
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FIG. 5: Chaplygin gas associated scalar field <j> evolution on the Bianchi type I brane universe as a function of the dimensionless 
time parameter r for 7 = 1/3, a = 0.5 and b = 0.005 (solid curve), b = 0.01 (dotted curve), b — 0.015 (dashed curve) and 
b = 0.02 (long dashed curve). We have assumed a vanishing effective cosmological constant on the brane, Ao = 0. The relation 
between r and the cosmic time t is t » 10 17 /i _1 x r s _1 , 0.5 < h < 1. 



1 

2? 



(52) 



flj = aio^ 1 ^ 6 ex P 



-1/4 



,1/2 



1,2,3, 



(53) 



12if 2 



(^) 



1/2' 



(54) 



a 2 = 



K 2 



(^) 



1/2 *' 



(55) 



9 = 5, 



(56) 



where we have assumed the initial conditions t = and V(0) = 0. 

In the opposite limit of small densities, the equation of state of the Chaplygin gas becomes p = —B/p a , corre- 
sponding to the limit 7 — > 0. In this case the equation describing the time dependence of the comoving volume V is 
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FIG. 6: The potential U(<j>) of the Chaplygin gas associated scalar field as a function of the scalar field <j> for 7 = 1/3, a — 0.5 
and b — 0.005 (solid curve), b = 0.01 (dotted curve), b = 0.015 (dashed curve) and b — 0.02 (long dashed curve). We have 
assumed a vanishing effective cosmological constant on the brane, Aq = 0. 



given by 



v= z Sv- 



B 



C 



B 



D 



C 



yi+a 



'~ 2,klB, 



2A 



-V 



B 



C 



yi+a 



(57) 



In the limit of large V the condition B » C/V 1+a is always fulfilled, and the general solution of Eq. ( 157(1 can be 
expressed as 



t-to = 



dV 



(58) 



where we have denoted 



A c = A 3fc2B 1 /(i+") l 



2A 



(59) 



In the same limit — > 00 we also have A C V 2 » 18U V 2 / 3 /kl\ + C x and therefore Eq. JSSJ can be immediately 
integrated to give 



= Vi exp [A c (f - t )] 



(60) 



where V\ is an arbitrary constant of integration. 

The time evolution of the physical and geometrical quantities in the limit of the large comoving volume is given by 



9 = K 



(61) 



a t = a i0 exp 



A,. 



(* - t ) 



exp 



cxp[-A c (t-i )] Lt = 1,2,3, 



(62) 



3K 2 

A = ^2 ex P i~ 2A c (* - *b) 



(63) 



if 2 

^ = 2V? ° XP [_2Ac ( * ~ to) 



(64) 



q = -l. 



(65) 
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In the limit of large times the mean anisotropy parameter tends to zero, A — > 0, and the Chaplygin gas filled 
Bianchi type I Universe ends in an isotropic state, corresponding to a de Sitter type Universe. This shows that the 
Chaplygin gas model offers a real possibility for replacing the effective cosmological constant on the brane. 
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